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Abstract 

We review limiting models for fracture in bundles of fibers, with statistically 
distributed thresholds for breakdown of individual fibers. During the break- 
down process, avalanches consisting of simultaneous rupture of several fibers 
occur, and the distribution D(A) of the magnitude A of such avalanches is 
the central characteristics in our analysis. For a bundle of parallel fibers two 
limiting models of load sharing are studied and contrasted: the global model 
in which the load carried by a bursting fiber is equally distributed among 
the surviving members, and the local model in which the nearest surviving 
neighbors take up the load. For the global model we investigate in particular 
the conditions on the threshold distribution which would lead to anomalous 
behavior, i.e. deviations from the asymptotics D(A) ~ A~2, known to be 
the generic behavior. For the local model no universal power-law asymptotics 
exists, but we show for a particular threshold distribution how the avalanche 
distribution can nevertheless be explicitly calculated in the large-bundle limit. 
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I. INTRODUCTION 



When a weak structural element in a material with stochastically distributed strengths 
fails, the increased load on the remaining elements may cause further ruptures, and thus 
induce a burst avalanche of a certain size A, i.e. one in which A elements fail simultaneously. 
When the load is further increased, new avalanches occur. The distribution of avalanche 
sizes, either at a fixed load, or the cumulative distribution from zero load until complete 
break-down of the material, depends on several factors, in particular the threshold strength 
distribution and the mechanism for load sharing between the elements. 

Due to the complex interplay of failures and redistributions of local stresses, few ana- 
lytical results are available in this field; computer simulations are commonly applied — see 
Herrmann and Roux [|l[ for a review. However, firm analytical results, albeit on simplified 
models, are important in order to develop a deeper understanding for universal properties 
and general trends. In the present article we therefore review and study burst events in 
models of fibrous materials that are sufficient simple to allow theoretical treatment. 

The models we consider are bundles of N parallel fibers, clamped at both ends, and 
stretched by a force F (Fig. [I]). The individual fibers in the bundle are assumed to have 
strength thresholds f\, i = 1,2, ... ,N, which are independent random variables with the 
same cumulative distribution function P(f) and corresponding density function p(f): 



Whenever a fiber experiences a force equal to or greater than its strength threshold, it breaks 
immediately and does not contribute to the strength of the bundle thereafter. The models 
differ, apart from differences in the threshold distribution, in how stress is redistributed on 
the surviving fibers when a fiber fails. 

A central quantity to be studied in the following is the expected number D(A,N) of 
bursts of size A when the fiber bundle is stretched until complete breakdown. 

The model of this kind with the longest history is one in which it is assumed that the 
fibers obey Hookean elasticity right up to the breaking point, and that the load distributes 
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(1) 



itself equally among the surviving fibers. The model with this democratic load redistribution 
is similar to mean-field models in statistical physics, and is called here the global model. 
For large iV Daniels || was able to determine the asymptotic distribution for the bundle 
strength, a result that been refined later . The distribution of burst avalanches was first 
studied by Hemmer and Hansen []5J . Their main result was that for a large class of threshold 
distributions P(f) the bursts were distributed according to an asymptotic power law, 



with a universal exponent 



lim — \-L ~ — , (2) 



* = §■ (3) 



In Sec. || we show that for special threshold distributions the power law (0) is not obeyed. 

The assumption of global load-sharing among surviving fibers is often unrealistic, and it 
is natural to consider models in which the extra stresses by a fiber rupture are taken up by 
the fibers in the immediate vicinity. The extreme version is to assume that only the nearest- 
neighbor surviving fibers take part in the load-sharing. In a one-dimensional geometry, as 
in Fig. [j], presicely two fibers, one on each side, share the extra stress. When the strength 
thresholds take only two values, the bundle strength distribution has been found analytically 
0. One interesting result is that the average bundle strength has a logarithmic size effect. 
The distribution of burst avalanches for such models with local load-sharing has not yet 
been determined, but simulations f7||| show that this model is not in the same universality 
class as the global model. The challenge to determine the burst distribution by other means 
than simulations remains, and that this is possible, at least in a special case (Sec. |TTTD , is 
one of the main results of the present article. 



II. THE GLOBAL MODEL 



In the global model the total force on a fiber bundle is distributed evenly on the surviving 
fibers. With a given load / per fiber, all fibers with threshold less than / have failed, while 
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fibers with thresholds above / survives. For large N the average number of fibers with 
thresholds exceeding the value / will be N[l — P(f)), corresponding to a average total force 

(F)(f)=Nf[l-P(f)]. (4) 

Let us for the moment assume that (F)(f) has a single maximum. This maximum corre- 
sponds then to the value f = f c for which d(F)/df vanishes. This gives 

l-P(fc)-fcP{fc)=0- (5) 

In Ref. the burst distribution was derived using the fiber elongation x as the indepen- 
dent variable, under the assumption that Hooke's law holds up to the threshold for breaking. 
Here, however, we formulate everything in terms of the force per fiber, /, and simplify the 
derivation by using directly the fact that the thresholds in a small interval of / are Poisson 
distributed. 



A. The burst distribution 

Let us consider a small force-per-fiber interval (/, / + df) in a range where the average 
force (F)(f) increases with /. For a large number N of fibers the expected number of 
surviving fibers is N[l — P(f)]- And the thresholds in the interval, of which there are 
Np(f)df, will be Poisson distributed. When N is arbitrary large, the burst sizes can be 
arbitrary large in any finite interval of /. 

Assume that an infinitesimal increase in the external force results in a break of a fiber with 
threshold /. Then the load that this fiber suffered, will be redistributed on the N[l — P{f)} 
remaining fibers; thus they experience a load increase 

,f = N\l-P(f)X (6) 
The average number of fibers that break as a result of this load increase is 

a = a(f) = Np(f) ■ 5f = J^fy (?) 

4 



For a burst of size A the increase in load per fiber will be a factor A larger than the 
quantity (||), and an average number a(/)A will break. The probability that precisely A — 1 
fibers break as a consequence of the first failure is given by a Poisson distribution with this 
average, i.e. it equals 

(A = l)! 6 • (8) 

This is not sufficient, however. We must ensure that the thresholds for these A — 1 fibers 
are not so high that the avalanche stops underway. This requires that at least n of the 
thresholds are in the interval (/, / + n5f), for 1 < n < A — 1. In other words, if we consider 
the A intervals (/, f + 5f),(f + 5f, f + 25 /), ...,(/ + (A - 1)5, f + A5f), we must find at 
most n — 1 thresholds in the n last intervals. There is the same a priori probability to find a 
threshold in any interval. The solution to this combinatorial problem is given in Appendix 
|X[ The resulting probability to find all intermediate thresholds weak enough equals 1/A. 
Combining this with (|8|), we have for the probability <f)(A, f) that the breaking of the first 
fiber results in a burst of size A: 

0(A,/) = ^a(/) A -V^. (9) 

This gives the probability of a burst of size A, as a consequence of a fiber burst due to an 
infinitesimal increase in the external load. However, we still have to ensure that the burst 
actually starts with the fiber in question and is not part of a larger avalanche starting with 
another, weaker, fiber. Let us determine the probability that this initial condition is 

fulfilled. 

For that purpose consider the d — 1 fibers with the largest thresholds below /. If there is 
no strength threshold in the interval (/ — 5f, /), at most one threshold value in the interval 
(/ — 25 f, /),... , at most d—1 values in the interval (/ — d5f, f), then fiber bundle can 
not at any of these previous /-values withstand the external load that forces the fiber with 
threshold / to break. The probability that there are precisely h fiber thresholds in the 
interval (/ — 5f d, f) equals 



h\ 

Dividing the interval into d subintervals each of length 5f, the probability ph,d that these 
conditions are fulfilled is exactly given by the solution of the combinatorial problem in 
Appendix A: ph,d = 1 — h/d. Summing over the possible values of h, we obtain the 
probability that the avalanche can not have started with the failure of a fiber with any of 
the d nearest-neighbor threshold values below /: 

WW = £ ^ e-d - = (1 - ^ + *>£ e-- (10) 

Finally we take the limit d — » oo, for which the last term vanishes. For a > 1 the sum must 
vanish since the left-hand side of ( |10D is non-negative, while the factor (1 — a) is negative. 
For a < 1, on the other hand, we find 

P b (f) = Jim A(/|d) = l-a, (11) 

a— »oo 

where a = The physical explanation of the different behavior for a > 1 and a < 1 

is straightforward: The maximum of the total force on the bundle occurs at f c for which 
a(/ c ) = 1, see Eqs. (^) and (0), so that a(/) > 1 corresponds to / values almost certainly 
involved in the final catastrophical burst. The region of interest for us is therefore when 
^ 1? where avalanches on a microscopic scale occur. This is accordance with what we 
found in the beginning of this section, viz. that the burst of a fiber with threshold / leads 
immediately to a average number a(f) of additional failures. 

Summing up, we obtain the probability that the fiber with threshold / is the first fiber 
in an avalanche of size A as the product 

Hf) = 0(A, f)P b (f) = ^ a(/) A - 1 e-^ A [l - o(/)], (12) 

where a(f) is given by Eq. (|7]), 

«(/) = J^J-y 
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Since the number of fibers with threshold in (/, f+Sf) is N p(f) df, the burst distribution 
is given by 

D(A) 1 rfc A A_1 rf° a 1 

N = NJo mp{f)df = ^rl a(/) 6 A [l-a(f)]p(f)df. (13) 

For large A the maximum contribution to the integral comes from the neighborhood 
of the upper integration limit, since a(f) e~ a ^' is maximal for a(f) = 1, i.e. for / = f c . 
Expansion around the saddle point and integration yields the asymptotic behavior 

D(A)/N oc A - !, (14) 

universal for those threshold distributions for which the assumption of a single maximum of 
(F)(f) is valid. 

Note that if the experiment had been stopped before complete breakdown, at a force per 
fiber f m < f c , the asymptotic behavior would have been exponential rather than a power 
law: 

D(A) /N OC A"* e - W/ m )-l-lna(/ m )]A_ ^ 

In the form 

J D(A)ocA-"e- A/Ao , with A Q oc (f c - f)~ v , (16) 

the breakdown process is similar to a critical phenomena with a critical point at total 
breakdown ||. The distribution follows a power law with index rj — | with a cutoff that 
diverges at total failure with an index v = ~. 

What happens when the average strength (F)(f) curve does not have a unique maxi- 
mum? If it has several parabolic maxima, and the absolute maximum does not come first 
(i.e. at the lowest / value), then there will be several avalanches of macroscopic size in 
the sense that a finite fraction of the iV fibers break simultaneously |p,0|| . The asymptotics 
(|T4]) is thereby unaffected, however. We turn next to threshold distributions that are more 
interesting because they lead to different asymptotics. 



B. Strong threshold distributions 



Rather than consider bundle strength functions (F)(f) with several parabolic maxima, 
we study now cases in which there is no such maximum. We are particularly interested in 
the asymptotics of the burst distributions. 

Model examples of such threshold distributions are 

I for / < /„ 

P(f) = (17) 
[l-[l + (/-/ )// r ]-« for />/„ 

Here a and fo are positive parameters, and f r is a reference quantity which we for simplicity 
put equal to unity in the following. This class of threshold distributions is rich enough to 
exhibit several qualitatively different avalanche distributions. 

The corresponding macroscopic bundle strength per fiber is, according to Eq. (|), 

(F)(f)_\f for f^fo (lg) 
N [ TTTTW for f > f° 
In Fig. |2| some threshold distributions p(f) and the corresponding macroscopic force curves 
(F)(f) are sketched. 

The distribution of avalanche sizes is given by Eq. (|13|). In the present case the function 
a(f) takes the form 

m fp(f) af / 1Q s 

T^Ptxj = TTT^To' (19) 

A simple special case is fo — 1, corresponding to 

p(f) = a r a ~ l for / > 1, 
since then the function ([19]) is independent of /: 

a(f) = a. 

This gives at once 



S 



D(A) _ 1 — a A A_1 r I- a . , 

N ~ ~a AT 



ae 



A-3 



(20) 



In other cases it is advantageous to change integration variable in Eq. ( |13D from / to a: 
D(A) A A - X 1 



J (a- o) a_1 (l - aK 1 (ae 1 "' 1 )' 1 da. (21) 



N e A A! a^-^l-Jo, 

3 

The asymptotics for large A, beyond the A~2 dependence of the prefactor, is determined 
by the A-dependent factor in the integrand. The maximum of ae 1_a is unity, obtained for 
a = 1, and the asymptotics depends crucially on whether a = 1 falls outside the range 
of integration, or inside (including the border). If the maximum falls inside the range of 
integration the D(A) oc A~5 dependence remains. A special case of this is a = 1, for which 
the maximum of the integrand is located at the integration limit and the macroscopic force 
has a "quadratic" maximum at infinity. 

Another special case is afo = 1 (and a < 1), for which again the standard asymptotics 
A~2 is valid. In this instance the macroscopic force has a quadratic minimum at / = /q 
(see Fig. |2| for a — 5), and critical behavior arises just as well from a minimum as from a 
maximum. 

In the remaining cases, in which a = 1 is not within the range of integration in Eq. ( pT| ) . 
the avalanche distribution is always a power law with an exponential cut-off, 

D(A) 



N 



A^ A A . (22) 



Here £ and A depend on the parameter values fo and a, however. This is easy to understand. 
Since 

<Hf) = «(! ~ fo) f2 oN 
df (I + /-/0) 2 ' 1 ) 

we see that a(f) is a monotonically decreasing function for fo > 1, so that the maximum of 
ae 1 ~ a is obtained at the lower limit / = / , where a = afo- The asymptotics 

D(A) oc A-l (a/ e 1-a/o ) A (24) 
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follows. 

This is true merely for afo < 1, however. For afo > 1 the macroscopic force (F)(f) 
decreases near / = /o so that a macroscopic burst takes place at a force /o per fiber, 
and stabilization is obtained at a larger force f\ (Fig. |2|). The subsequent bursts have an 
asymptotics 

D(A)ocA-|(a(/ 1 )e 1 ^)) A , (25) 

determined by the neighborhood of / = j\. 

For /o < 1, the maximum of ae l ~ a is obtained at / = oo, leading to the asymptotics 

D(A) oc A-i- a (ae 1 " a ) A , (26) 

reflecting the power-law behavior of the integrand at infinity. 

The results are summarized in Table |. Note that the /o = 1 result (p0|) cannot be 
obtained by putting /o = 1 in Eq. ( p4| ) since in fl2T|) the order of the limits A — > oo and 
/o — > 1 is crucial. 



III. THE LOCAL MODEL 

The assumption of global loadsharing among surviving fibers is often unrealistic, since 
fibers in the neighborhood of the failed fiber are expected to take most of the load increase. 
The extreme form for local load redistribution is that all extra stresses caused by a fiber 
failure are taken up by the nearest-neighbor surviving fibers. 

The simplest geometry is one-dimensional so that the N fibers are ordered linearly, 
without or with periodic boundary conditions (Fig. [1]). In this case precisely two fibers, one 
on each side, take up, and divide equally, the extra stress. At a total force F tot on the bundle 
the force on a fiber surrounded by ni previously failed fibers on the left-hand side, and n r 
on the right-hand side, is then 

I !j±(l + \(n l + n r ))=x{2 + n l + n T ). (27) 

10 



Here 

one-half the force-per-fiber, is a convenient variable to use as the driving force parameter. 
This model has been discussed previously [llj for a different purpose. Preliminary studies 



|7||| of the avalanche distribution for some threshold strength distributions have not yielded 
analytical results but simulation results that show convincingly that the local model is not 
in the same univerality class as the global model. 

In order to obtain explicit results we assume for the fiber strengths the simplest possible 
case, a uniform threshold distribution. In units of the maximum threshold: 

I / for < / < 1 
P(f) = ~ (29) 

I 1 for / > 1. 

Avalanches in the local and the global models have different characters. In the local 
model an avalanche unroll with one failure acting as the seed. If many neighboring fibers 
have failed, the load on the fibers on each side is high, and if they burst the load on the 
new neighbors will be even higher, etc. In this way a weak region in the bundle may be 
responsible for the failure of the whole bundle. For a large number N of fibers the probability 
of a weak region somewhere is high, and this explains in a qualitative way that the maximum 
load the bundle are able to carry does not increase proportional to N, but slower than linear. 

The load distribution rule (|27|) implies that an avalanche of size A does necessarily lead 
to a complete breakdown of the whole bundle if the external force is too high, i.e., if x 
exceeds a critical value x max . Since here a fiber can at most take a load of unity, we have 

1 

Xmax — ^ — 2 • 

The strategy of the derivation is to first establish a set of recursion relations between 
quantities that give probabilities of certain configurations at fixed external force, i.e., at 
fixed x. Afterwards (Sec. [Ill rj| ) we connect this with the size distribution of avalanches for 
all x up to the critical value x max . 
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A. Recursion relations 



We will use the terminology that the magnitude of an avalanche is the number of failing 
fibers in the avalanche, and the length of an avalanche is the number of fibers between the 
nearest surviving fibers on each side of the avalanche. The length can be larger than the 
magnitude since it may include fibers that have failed in previous avalanches. 

We define S(l; x), the gap probability, to be the probability (at given force parameter x) 
that in a selected region of I consecutive fibers all fibers have failed, assuming the two fibers 
on each side to be intact. We let S(0; x) = 1 by definition. 

Another central quantity is the probability density p(l,a;x). We define it by selecting 
a region of I consequtive fibers, and let p(l, a; x) dx be the probability that a force increase 
from x to x + dx leads to an avalanche of this length I and of magnitude a. 

The state at force parameter x that all L fibers have failed must have appeared for some 
force parameter in the range (0, x), and by a burst of some magnitude a in the range (1, L). 



Let us now obtain expressions for the probability density p(L, a;x), first for the special 
case that the magnitude a is unity. Just one fiber fails in this burst, and in an avalanche of 
length L therefore L — 1 of the neighboring fibers must already have failed. By fl2~7|) the force 
on the fiber just before it fails is (L + l)x, and for the uniform distribution the probability 
that it fails due to a force parameter increase dx is just (L + l)dx. The probability of the 
burst of magnitude 1 to occur when x —>■ x + dx is this probability of failure of the single 
remaining fiber, (L + l)dx, times the probability that the L — 1 neighbors have already 
failed. The latter is given by the appropriate gap probabilities. Since the position of the 
failing fiber is arbitrary in the interval we have 



Thus 




a=l 



(31) 



L-l 




E 




(32) 
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We next consider expressions for the probability p(L, a; x) with an internal avalanche 
of magnitude a larger than unity. For that purpose we introduce two new quantities: Let 
Pi(L + l,a;x) dx be the probability that a fiber fails because a force parameter increase 
x — > x + dx starts, on its right-hand side, an avalanche of magnitude a (not counting the 
ultimate fiber on the left-hand side) and of length L. Similarly p r (L + l,a;x) dx is the 
probability that a fiber fails because the force parameter increase x — > x + dx starts, on its 
left-hand side, an avalanche of magnitude a and of length L. 

Consider the event described by p(L, a; x), and let the last of the a fibers that fail be fiber 
T. The force distribution mechanism in the local model implies that T is either the leftmost 
or the rightmost of the a fibers. The first possibility implies that the increase x — > x + dx 
induces the first failure to the right of JF, which starts an avalanche of magnitude a — 1 and 
length i, say, to the right of T. Here a — 1 < i < L — 1, of course, and T must have L — i — 1 
previously failed fibers on its left-hand side. 

Including all possibilities we have 

L-l 

p(L,a;x) = ^2 [S(L — i — l;x) pi(i + l,a— l;x)+ 

i=a—l 

p r (i + 1, a — 1; x) S(L — i — 1; x)] , (33) 

where the second term represents events in which the first failure occurs to the left of T . 

On the other hand we want to express pi(l+l, a; x) andp r (/+l, a; x) in terms of previously 
defined quantities. For magnitude a = 1 this is relatively simple. Let the single fiber (call 
it Q) that starts the process have n already failed fibers to the right and I — n — 1 fibers 
to the left, with < n < I — 1. The probability that fiber Q fails under the load increase 
x — > x + dx is (I + l)dx for the uniform threshold distribution. Since the failure of fiber Q 
causes a load increase (n+ l)x on the left fiber, the probability of its failure is (n + l)x. We 
have, when all possible positions of Q are taken into account, 

i-i 

pi(l + 1, 1; x) dx = S ( l ~ n ~ !; x ) ( l + l ) dx S ( n 'i x )( n + l ) x ( 34 ) 

n=0 

Similarly, 

13 



l-l 

p r (l + l,l;i) = ^S(l-n- 1; x)(Z + l)S(n; a:)(Z - n)x = p,(Z + 1, 1; x). (35) 

n=0 

The corresponding expressions for p r (L + 1, a; x) and p t (L + l,a;x), with a larger than 
unity, are more complicated. The internal avalanche started by x ^ x + dx proceeds so that 
the final failure is either the leftmost or the rightmost of the a fibers, or both. If both go 
simultaneously, we make the arbitrary definition that in such a case the right-hand neighbor 
fails first. This secures a unique sequential ordering of failures. 

Consider first p r (L + 1, a + 1; x), and denote the last surviving fiber on the right-hand 
side as T . Assume first that the rightmost of the a internal fibers fails last, and let this fiber 
have i fibers on its left-hand side and L — i — 1 failed fibers on the right-hand side. The 
probability that this right-hand side fiber fails under x — > x + dx with an internal avalanche 
of magnitude a and length i is just p r (i + 1, a; x) dx, and the probability of finding L — i — 1 
failed fibers on the right-hand side is given by the gap probability S(L — i — l;x). After the 
rightmost internal fiber has failed the load increase on the T is [i + l)x, which also equals 
its probability of failure. The other alternative is that the leftmost of the a internal fibers 
fails last, with, say, i fibers on its right-hand side, and L — i — 1 failed fibers on its left. Then 
the extra load increase on J 7 , and hence its probability of failure, is (L — i)x. Including all 
possible positions i we end up with 

L-l 

p r (L + 1, a + 1; x) — S(L — i — 1; x) \p r (i + 1, a; x)(i + l)x+ 

Pi(i + 1, a; x)(L — i)x] . (36) 

By a similar argument the corresponding expression for pi(L + 1, a + 1; x) is built up. 
However, when in this case the rightmost of the internal fibers fails last we must add the 
probability that the leftmost and the rightmost fibers are simultaneously overburdened. 
Letting p 2 (i + 2, a; x) dx be the probability that an avalanche of length i and size a makes 
both neighbor fibers fail, we have 

L-l 

Pi(L + 1, a + 1; x) — S(L — i — 1; x) \pi(i + 1, a; x)(i + l)x+ 

i=0 

p r (i + l,a;x)(L - i)x + p 2 (i + 2,a;x)] . (37) 
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Finally, the recursion relations for p 2 close the set. One sees easily that when the failure 
of the two end fibers is caused by a single internal fiber burst due to the force increase from 
x to x + dx, we have 

L-l 

p 2 {L + 2, 1; x) = ^ S(L - i - !; X )( L + l ) S & + l)x(L - i)x. (38) 

i=0 

Here (L + l)dx is the probability that the single fiber fails, (i + l)x is the probability that 
the (i + 1) new failures on the right makes the left-hand-side fiber break, while (L — i)x 
is the probability that the [L — %) new failures on the left makes the right-hand-side fiber 
break. 

When the failure of the two end fibers are caused by an internal avalanche involving 
a > 1 fibers we may argue along the same lines as for pi, with the result 

L-l 

p 2 {L + 2, a + 1; x) — S(L — i — 1; x) [pi(i + 1, a; x)(i + l)x(L — i)x 

i=0 

+ (p r (i + 1, a; x)(L - i)x + p 2 (i + 2, a; x)) (i + l)x] . (39) 

We can simplify the set of equations somewhat by introducing the sum p s = pi+p r , with 
the result 



S(L; x 
p(L, l;x 
p(L, a + 1; x 
p s (L + 1, l;x 
p s {L + 1, a + 1; x 

p 2 (L + 2,l;x 

p 2 (L + 2, a + l;x 



L i-x 

a=l J ° 

L-l 

J2S(L-i- l;x)(L + l)S(i;x) 

i=0 
L-l 

^ S(L - i - 1; x)p s (i + 1, a; x) 

i=0 
L-l 

x)(L + 1)5(1; x)(L + l)x 

i=0 
L-l 

^ 5(L - i - 1; x) \p s (i + 1, a; x) 

i=0 

+p 2 (i + 2,a;x)] 

L-l 

^5(L-i- l;x)(L + l)5(i;x) 

i=0 

(i + l)x(L — i)x 

L-l 

5(1/ — i — 1; x)(i + l)x [p s (i + 1, a; x)(L — i)x 

i=0 

+p 2 (i + 2,a;x)] 



(40) 
(41) 
(42) 
(43) 

(44) 
(45) 
(46) 
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Starting with S(0; x) = 1, one easily proves by induction the following x-dependence of 
all quantities involved: 

S{L; x) = S(L)x L (47) 
p(L,a;x) = p(L,a)x L ~ 1 (48) 
p s (L + l,a;x) =p s (L + l,a)x L (49) 
p A (L + 2,a;x) =p A (L + 2,a)x L+1 (50) 

In ^-independent form the recursion relations then take the form (with a > 0) 

S{L) = L- l j^p{L,a) (51) 

a=l 

p(L, 1) = J2S(L- i - 1)(L + l)S(i) (52) 

i=0 
L-l 

p(L, a + 1) = S ( L ~ * - !)P-(i + !> a ) (53) 

i=0 
L-l 

p s (L + l,l) = y £S(L-i- 1)(L + l)5(z)(L + 1) (54) 

i=0 
L-l 

p s (L + 1, a + 1) = J2 S i L ~ i ~ !) \Ps(i + 1, a) + P 2 (i + 2, a)] (55) 

i=0 
L-l 

p 2 (L + 2, 1) = S{L - i - 1)(L + + 1)(L - i) (56) 

8=0 

L-l 

j9 2 (L + 2,a+ 1) = ^ S(L-i-l)(i + l) \p s (i + l,a)(L-i)+ 

i=0 

p 2 (z + 2,a)] (57) 

Let us finally note that the feature of the uniform distribution that makes the derivation 
simpler than for other distributions is that the probability for failure of a fiber is given by 
the load increase, independent of the actual load level. 

We can now calculate recursively S(L; x) and p(L, a; x) for integer L and a. By d47| ) -(|50 D 
the x dependence is trivial. 
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B. The asymptotic burst distribution 



In order to use the quantitative information obtained above we must first determine the 
survival probability P s (N,x) that a fiber bundle is able to tolerate a force per fiber equal 
to 2x. Noting that in this model avalanches are local phenomena, and that two failed fibers 
are only correlated when all fibers inbetween have failed, the survival probability P s (N,x) 
is expected to depend exponentially on the length N for large N, so that 

lim N' 1 In P S (N, x) = t(x) (58) 

AT— >oo 

is finite. The exponential form of the survival probability is discussed and confirmed in other 
studies 0-0. 



We assume periodic boundary conditions, and number the fibers from an arbitrary start- 
ing point. We define Pf(n, L; x) to be the probability, at force parameter x, that among the 
n first fibers there is no fatal burst, and that the last L fibers of these have all failed. Fiber 
number n + 1 is assumed to hold. We will now establish a recursion relation between the 
P f (n,L;x). 

Consider a region of n + 1 + L fibers in which no fatal burst has occurred, and where 
the last L fibers have failed. The probability of this configuration is Pf(n + 1 + L, L; x). In 
the region to the left of fiber number n + 1 let the length of the region of broken fibers that 
contain fiber number n be i, where i may take all values between zero (if fiber number n is 
intact) and M(x) = [or 1 - 2]. The region to the right of fiber number n + 1 has L broken 
fibers, and the probability of this is S(L; x). This gives the recursion relation 

M(x)-L 

P f (n + l + L,L;x) = J2 Pf{n,i;x)[l - {i + L + 2)x]S{L;x). (59) 

i=0 

The last factor is the probability that fiber number n+1, which has i + L failed neighbors, 
holds. 

Insertion of the product form Pf(n,L;x) ~ t(x) n Pf(L; x) into (|59|) yields the following 
equations for the Pf(i; x): 
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MO) 

P f (L; x)- £ [1 - (i + L + 2)x]S(L] x)t(x)^ 1 P f {i; x). (60) 

It is consistent to let P/(0;x) = 1. Since L may take the values 0, 1, . . . ,M(x), ( |60"D is a 
set of M(x) + 1 homogeneous equations for the M(x) + 1 quantities Pf(i; x). The system 
determinant of the equation set must vanish, and this determines t(x) for a given force 
parameter x. With Pf(0;x) = 1 all quantities can then determined. The practical solution 
procedure is by iteration. 

From the definitions of Pf(n, L; x) and S(L; x) it follows that the ratio 

Pf(n,L\x) 
S{L-x) 

is the probabilility, at force parameter x, that among the first n fibers there is no fatal burst, 
given that there are L\ failed fibers on the right-hand side. Then 



P f {n,L x ;x) 



p(Li, A; x) dx (61) 



S(L;x) 

is the probability that an increase of the force parameter from x to x + dx starts an avalanche 
of size A and length L, so that afterwards there is no fatal burst among the n fibers on the 
left-hand side. 

Finally we want to determine the probability for a burst of size A in a system of N fibers 
in a ring configuration (Fig. [l]). On the left of a selected fiber f we consider a region of n 
fibers, and on the right a region of iV — n — 1 fibers. The probability that the force parameter 
increase x — * x + dx induces a burst of size A and length L\ to the left of f that holds is 
given by (|6l|) . Here A < L\ < n, of course. On the right-hand side of f a number L2 fibers 
adjacent to f may have failed. (Here L 2 is less than the remaining number iV — n — 1 of 
fibers.) The probability of such a configuration (with no fatal burst) is Pf{N — n — 1, L2; x). 
We must also take into account that the fiber f itself, with Li + L 2 failed neighboring fibers, 
must hold, the probability of which is [1 — (Lx + L 2 + 2)x]. 

When we take this together, sum over the possible values of L\, L 2 and n, and integrate 
over x, we obtain 



l/(A+2) N M{x)M{x) _ Ll 

d(a)= / E E E siL- f ^^ 

J Q n=lL 1= A L 2 =0 0[L,,X) 

P f (N-n-l,L 2 ;x)[l-(L 1 + L 2 + 2)x] dx. (62) 

Using the product property Pf(n,L;x) ~ t(x) N Pf(L; x) the sum over n simply yields a 
factor N, and we find 

l/fA+2) 

[1 - (L x + L 2 + 2)z]dx. (63) 

This may now be evaluated. The results for a bundle of N = 20 000 are shown in Table 
HI together with simulation results for 4 000 000 bundles, each having 20 000 fibers. 

The agreement between the simulation data and the theoretical data is, as we see, ex- 
tremely satisfactory. 

An analysis of the burst distribution obtained for this local model shows that the distri- 
bution does not follow a power law except for small values of A (Fig. |3|). If one nevertheless 
does a linear regression analysis on this part of the data set, the effective power would be 
of the order 5, considerable larger than the "mean-field" value | for the global model [|7|||. 

C. Size-dependent bundle strength 

Let us now attempt to find an simple estimate for the maximal force per fiber that the 
fiber bundle can tolerate. In order to do that we assume that the fatal burst occurs in a 
region where no fibers have previously failed so that the burst has the same magnitude and 
length. We know that a single burst of length A = x~ x — 2 is fatal, Eq. (|30D, so our criterion 
is simply 

Z^aT 1 - 2) = 1. (64) 

If we take into account that the two fibers adjacent to the burst should hold, and ignore 
the rest of the bundle, the gap distribution would be 
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l/(A+2) 



N~ 1 D(A) w /" [l-(2 + A)ifp(A,Aix)& 



2p(A,A) 



A(A + 1)(A + 2) A + 1 ' 



(65) 



With the abbreviation 



Rt 



P(A,A) 
(A-l)!' 



we have 



D(A)/N 



2(A + 2)! 



A 2 (A + 1) 2 (A + 2) A + 2 

using Stirling's formula. 

Taking logarithms we have 

lni? A 



8tt(A + 2) 
A 2 (A + l) 2 



-A-2 



Ra, 



(66) 



lnD(A) -]nN= -(A + 2) 



1 + 



A + 2 



+ o 



'lnA 



-(A + 2), 



(67) 



using the result ( p9|) of Appendix |B| for when A is large. 
The failure criterion (1641) then takes the form 



In 



x 



Since x = F/2N we have the following estimate for the maximum force F that the fiber 
bundle can tolerate before complete failure: 

Thus the maximum load that the fiber bundle can carry does not increase proportional to 
the number of fibers, but slower. This is to be expected since the probability of finding 
somewhere a stretch of weak fibers that start a fatal avalanche increases when the number 
of fibers increases. 



The N/ In iV dependence agrees with a previous estimate by Zhang and Ding |L5[ and is 



seen also in the model with thresholds zero or unity flq,|12|| 
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IV. CONCLUDING REMARKS 



We have in this article discussed burst distributions in fiber bundles with two different 
mechanisms for load distribution when fibers rupture, viz. global or extremely local load 
redistributions. 

The main results are the following: 

(i) For the global model the burst distribution follows a universal power law A~5. 

(ii) Deviations from this power-law dependence may, however, occur for exceptional distri- 
butions of fiber strengths. 

(iii) For the local model and for a uniform distribution of fiber thresholds we show that 
it is possible, although complicated, to carry through an theoretical analysis of the burst 
distribution. 

(iv) A simulation study for a bundle of 20 000 fibers confirms convincingly the theoretical 
results. 

(v) For the local model the burst distribution falls off with increasing burst size much faster 
than for the global model, and does not follow a power law. 

(vi) The expected maximum load that a bundle with global redistribution mechanism can 
tolerate increases proportional to the number N of fibers, and proportional to N/ In N for 
the local redistribution mechanism. 



APPENDIX A: 



The combinatorial problem in Sec. |II A| can be formulated more generally as follows: Let 
Ph t n be the probability that by distributing h nonidentical particles among n numbered boxes, 
box number 1 will contain no particles, box number 2 will contain at most 1 particle, and in 
general box number i will contain at most i — 1 particles. 

Since the probability that there are h — k particles in box number n is equal to 

(h\ n\ h ~ k fn-l^ k 



k J \nJ V n 
21 



we must have 



We now prove by induction that 

Ph,n = 1 - -• (A2) 

n 

Assume that this holds for ph, n -i, ah h. Insertion into the right-hand side of (|A1|) gives 

'h\ n\ h - k /n-V ' 



Ph,n = 



k J \nj \ n J \ n — 1 J 

in accordance with (|A2|). Since is valid for n = 2, the induction is complete. 
For the application in the text, 

1 

Pn—l,n 

n 

is needed. 



APPENDIX B: 

In Sec. |III C| an estimate for p(n, n, ) was needed. We base it on the recursion relations 
(|53|), ( j55|) and (|57D for L = n + 1, a = n: 

p(n + 1, n + 1) = p s (n + 1, n) (Bl) 
p s {n + 2, n + 1) = {n + 2)p s (n + 1, n) + p 2 {n + 2, n) (B2) 
p 2 {n + 3, n + 1) = (n + l)p s (n + 1, n) + (n + l> 2 (n + 2, n) (B3) 

We have used that p s (n + 1, a) og + 2, a) vanish for a > n. 
It is easy to eliminate p s by ([Bl]) and p 2 by ( p2|) , with the result 

p[n + 1, n + 1) = 2np(n, n) — (n — l) 2 p(n — 1, n — 1). (B4) 
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This is a three-term recursion starting off with p(l, 1) = 2 and p(2, 2) = 2p(l, 1) by ( |B4|) . 
With 



Rn 



p(n, n) 
(n - 1)1 



the recursion takes the form 



(B5) 



Rn+l — 2R n — (1 _ -)Rn-l- 



(B6) 



Introducing the generating function 



G(z) = £ R n z n 



n=l 



the recursion (p6|) may be transformed to the differential equation 





dz 



G(z)(l-z 2 )/z =G(z) 



with solution 



G(z) 



2z 
1 - z 



,*/(!-*) 



(B7) 



(B8) 



Thus the radius of convergence of the power series ( |B7D is unity and therefore 



lim Ri /n = 1. 



(B9) 



In fact R n oc n le 2 ^™ for large n [|T| 
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FIGURES 

FIG. 1. A fiber bundle with periodic boundary conditions. 

FIG. 2. The threshold distribution density p(f) and the macroscopic bundle strength (F)(f) 
for the distribution (17), with fo = 2f r , and for a = | (upper curve), \ (middle curve), and 
| (lower curve). The broken part of the a = 2/3-curve is unstable and the macroscopic bundle 
strength will follow the solid line. 

FIG. 3. Burst distribution in local model as found numerically for 4 000 000 samples with 
N = 20000 fibers (+), and calculated from Eqs. (68) and (69) (o). The straight line shows the 
power law A -5 and the broken curve the function exp(— A/Aq) with Ao = 1.1. Note the small 
value of Aq. 
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TABLES 



Parameters Asymptotics 

< f < 1, a < 1 A-§- a (ae 1_a ) A 

0</ <l,a = l A _ i 

/o = l,a<l A-i(ae 1 - a ) A 

l<fo<a~ 1 A"f (a/ e 1 - a/o ) A 

1< f = a- 1 A-f 

1< a- 1 < f Q A-f e- A / A « 
TABLE I. Asymptotic behavior of the burst distribution for strong threshold distributions in 
the global model. 
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A Simulation Calculation 

1 8 327 378 752 8 327 331 808 

2 491 305 573 491 331 178 

3 72 126 803 72 114 644 

4 17 179 080 17 180 414 

5 5 590 887 5 591 243 

6 2 243 916 2 243 012 

7 1 030 833 1 031 678 

8 515 309 515 310 

9 268 589 268 139 

10 140 911 140 751 

11 72 251 72 701 

12 36 525 36 277 

13 17 523 17 285 

14 8 015 7 835 

15 3 352 3 392 

16 1 442 1 418 

17 559 579 

18 223 233 

19 90 93.8 

20 40 37.5 

21 18 15.0 

22 10 6.0 

23 1 2.4 

24 2 1.0 

25 0.4 
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TABLE II. The burst distribution D(A) for the local model with a bundle of N = 20000 
fibers. The simulation results are based on 4 000 000 samples. The calculated values are based on 
Eqs. (68) and (69), and have been multiplied by 4 000 000. 
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